Phase diagram of ultracold atoms in optical lattices: Comparative study of slave 
fermion and slave boson approaches to the Bose-Hubbard model 
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We perform a comparative study of the finite temperature behavior of ultracold Bose atoms in 
optical lattices by the slave fermion and the slave boson approaches to the Bose Hubbard model. 
The phase diagram of the system is presented. Although both approaches are equivalent without 
approximations, the mean field theory based on the slave fermion technique is quantitatively more 
appropriate. Conceptually, the slave fermion approach automatically excludes the double occupancy 
of two identical fermions on the same lattice site. By comparing to known results in limiting cases, 
we find the slave fermion approach better than the slave boson approach. For example, in the non- 
interacting limit, the critical temperature of the sup erfiuid- normal liquid transition calculated by 
the slave fermion approach is closer to the well-known ideal Bose gas result. At zero-temperature 
limit of the critical interaction strength from the slave fermion approach is also closer to that from 
the direct calculation using a zero-temperature mean field theory. 
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I. INTRODUCTION 

Strongly correlated systems are of longstanding inter- 
est in studies of condensed matter physics. Ultra-cold 
atoms in optical lattices offer new opportunities to study 
strongly correlated phenomena in a highly controllable 
environment ^ ,2;i 4J ■ A quantum phase transition, the 
superfluid/Mott-insulator transition, was demonstrated 
using *^Rb atoms in three- and one-dimensional lat- 
tices 4|. Strongly correlated phenomena for boson sys- 
tems may be studied theoretically by the Bose-Hubbard 
model p. Experimental feasibility was established by 
microscopic calculations of the model parameters for cold 
boson atoms in optical lattices A review of recent 
works on the superfluid-insulator quantum phase transi- 
tion at zero temperature is given in ref. 0- 

Strictly speaking, a quantum phase transition can not 
be observed at any finite temperature. The experimen- 
tal data give only a signal that the system is towards a 
quantum phase transition if the temperature is extrap- 
olated to zero. What the experiments really observed 
was a transition from the superfluid to the normal liquid 
whose compressibility is very close to zero and the system 
is practically a Mott insulator. Such a 'classical' phase 
transition has been investigated recently by Dickerscheid 
et al Phase diagrams for a given atom density were 
calculated in the temperature-interaction plane and the 
chemical potential-interaction plane: For a commensu- 
rate optical lattice, there are only the superfluid and the 
Mott insulator phases at zero temperature. At finite tem- 
peratures, starting from the superfluid phase, there is a 
superfluid/normal liquid phase transition while the Mott 
insulator phase crossovers to the normal liquid. 

In order to extend the ordinary mean flcld approach for 
the Bose Hubbard model _9J to include the flnite tempera- 
ture effects, the slave boson technique ^3 used. The 
slave particle technique has been widely applied in deal- 



ing with the strongly correlated electron systems ■ In 
principle, the slave boson and slave fermion approaches 
are equivalent. However, in practical calculations, ap- 
proximations still have to be used. It was well-known 
that in the t-J model, the same mean field approximation 
using the slave boson or slave fermion leads to very differ- 
ent phase diagrams ■ K was known that the slave bo- 
son mean field approximation can qualitatively describe 
the phase diagram of the cuprates at finite doping. How- 
ever, due to the Bose condensation of the holons, the 
slave boson mean field approximation does not produce 
correctly the ferromagnetic Mott insulator phase. One of 
the purpose of this work is to examine if both slave parti- 
cle approaches to the Bose Hubbard model give the same 
physical results under the mean field approximation. 

In these approaches, there is a constraint that each 
site can be occupied by only one slave particle. With 
this exact constraint, the model is very hard to solve. A 
standard approximation is to relax the constraint on each 
site to the requirement that the average slave particle per 
site over the lattice be equal to 1. While both approaches 
give the same qualitative phase diagram, we shall see that 
the quantitative behaviors derived from the slave fermion 
approach are more accurate. The advantage of the slave 
fermions is that the Fermi statistics automatically ex- 
cludes two same type slave fermions from occupying the 
same site even when the constraint is relaxed. We shall 
see that the configurations with the multi-occupations 
of different types of slave fermion are far away from the 
mean field state we consider. Thus, these configurations 
will not significantly influence our results. However, the 
statistics of the slave particle affect the result remarkably. 
For repulsive interactions, there are two unsatisfactory 
features arising from the finite temperature mean field 
theory in the slave boson approacheQ- One of them is 
that the critical on-site repulsive Uc ~ 5.83 from a direct 
zero temperature mean field calculation and differs from 
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the critical ?7c ~ 6 in T — > from the finite temperature 
mean field. This difference is much smaller with the slave 
fermion approach. The other was that there is a maxi- 
mum Tc at U =/= in the U-Tc curve which is obviously 
unphysical. Wc find that these deficiencies are corrected 
in the slave fermion approach. Furthermore, for U ^ 0, 
the critical temperature of the superfluid-normal liquid 
phase transition for the ideal Bose gas was well-known. 
We find that this critical temperature calculated by the 
slave fermion approach is much closer to its exact value 
than that by the slave boson approach. 

This paper is organized as follows. In Sec. II, we give 
an overview of our slave particle approaches. In Sec. Ill, 
the perturbation theory is introduced. In Sec. IV, we 
give our main results according to the mean field theory. 
Section V is our conclusion. 
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II. SLAVE PARTICLE APPROACH 

A boson operator on site i may be expressed by the 
occupation state |a), i.e.. 



a| = ^ Va + l\a + l)ii{a\. 



(1) 



The slave particles are the auxiliary particles, which 
are obtained by mapping the occupation state \a)i 
aa,i- According to this mapping, the original boson cre- 
ating operator at on the lattice site i may be decomposed 

as 



(2) 



where Ua^i may be either the (slave) boson operator ba.i 
with [ba,i,bj^j] = SapSij or the (slave) fermion operator 

Ca.i with {ca,i, Cp j} — SapSij. As the auxiliary particles, 
they have to obey the constraint 



(3) 



on each site, which corresponds to the completeness of 
the states - J^a ~ ^ ^^'^ ^^'^ original Bose commu- 

tation relation: [0^,0]] = Sij. 

The Bose Hubbard Hamiltonian we will focus on reads 

H = -t^aTaj - /i^n,; + j'^n,{n, - 1), (4) 

where the symbol {ij) denotes the sum over all nearest 
neighbor sites, is the chemical potential. The hopping 
amplitude t and the on-site interaction U are defined by 

r 

t= drW*(r)\ + V(r)]W(r + a), 

J 2m 

U ^ g j dv\W{v)\\ (5) 



FIG. 1: The dependence of the Bose Hubbard parameter 
aaU/Qt for the three dimensional optical lattice as a function 
of the strength of the lattice potential Vq. The scattering 
length as is in units of a nanometer. Vb is in units of the 
recoil energy 



where g 



with as the s-wave scattering length of 



the atoms and m the mass of atom; W{v) is the Wan- 
nier function corresponding to the lowest Bloch band and 
V{r) is the periodic optical lattice potential. In three di- 
mensions, taking the periodic potential with a form 

V{v) = Va{sm^{kx) + iiin^{ky) + sm^{kz)), (6) 

we can calculate U and t from the band theory 0. In 
Fig. 1, we display the relation between U/Qt and Vq. 

In the slave particle language, the one-component Bose 
Hubbard Hamiltonian in a d-dimensional cubic lattice 
with L sites reads 



- A.^5]a< + ^^^a(a~l)<, (7) 

i a i a 

With the constraint ((^J, the standard path integral leads 
to the partition function of the system 



U 



—a{a - 1) - iXi]aa,z + 



(8) 



(ij) a(3 

where we have taken h = ks = 1- dai is equal to the 
complex conjugate &* ^ of the bosonic field ba^i for the 
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slave boson while being the Grassmann conjugate Cai of 
the fermion field c^i for the slave fermion. The integrals 
over the Lagrange multiplier field Ai(r) come from the 
constraint 



:(^<-l)dr 



(9) 



In the sense of the (5-function, the Xi fields have to be real 
to ensure the constraint is correctly taken into account. 

To decouple the four slave particle term in the Hamil- 
tonian, we introduce a Hubbard- Stratonovich field 
which is a bosonic field and may be identified as the order 
parameter of the superfluid. The integral 

I cxp - 1 drtY^ 

($* - ^ Va + laa+i,iaQ,i) 

a 



(10) 



is obviously a constant. The partition function can be 
written as 



U 



— a{a - 1) - i\i]aa^i + i ^ Ai 

i 



(11) 



(12) 



This effective action is the starting point of the slave par- 
ticle approach. The slave boson and the slave fermion 
representations are only reexpression of the original Bose 
Hubbard model. Both of them are equivalent to the orig- 
inal model before any approximation. So far, all formal 
transformations we made are rigorous. 



III. PERTURBATION THEORY 

The physical meaning of the $ field may be seen by 
using its equation of motion: 



(13) 



This means that $i indeed serves as an order parameter 
field. Near the Mott transition, this order parameter is 



small and one can use perturbation theory to solve the 
system described by the action ((T^ . The difficulty is 
that there is no way to exactly solve the problem if the 
A field varies from site to site. A widely-used approx- 
imation is to relax the constraint (3) by replacing the 
local constraint Lagrange multiplier Ai(r) by an imagi- 
nary time- and site-independent field A. That is, relax- 
ing the condition of exactly one slave particle per site 
to one with an average of one particle per site. It im- 
plies that multi-occupation of the slave particles on the 
same site is allowed. For slave bosons, this relaxation al- 
lows the same type of the boson to multi-occupy a single 
site. However, for the same type of slave fermions, multi- 
occupation of the same site is automatically forbidden by 
the Pauli principle. The value of A will be variationally 
determined. 

To do the perturbation calculation, it is convenient to 
make a Fourier transformation for the fields Ai ^ Ui, $i 
and A, : 



Ai. 



1 



/LP 



:-\ — lUJnT 



(14) 



k,n 



where the Matsubara frequencies LUn = 2TmT for bosonic 
fields and a;„ = (2n + 1)ttT for fermionic fields. The 
approximation of the site-independent of A^ implies that 
all Xk,n = Ao,o = Xy/Lp. 

After all these preparations and then some algebra, we 
arrive the effective action to second order of the order 
parameter field, 

^£,cff[$*,*] =/3f^o-^$k,„G"'(k,zc^„)$k,„, (15) 

k,n 

where $k.n is the Fourier component of (E>i and the zero 
order thermodynamic potential is given by 



(16) 



where /3 — 1/T and eo(a) = —iX — afi + a{a — l)U/2 
and — (-I-) sign corresponds to the slave boson ( slave 
fermion) approximation. The Green's function is defined 
by 



[I — all 



.(17) 



where the dispersion ek = 2i ^ • cos(fcio) with z the near- 
est neighbor partition and a the lattice spacing vector of 
a c?-dimensional cubic lattice. The bosonic Matsubara 
frequency a;„ = . The slave particle occupation num- 
ber is given by 



(18) 



exp{/3[-iA - a/i + a{a - l)U/2]} =F 1 ' 



corresponding to the slave boson and slave fermion, re- 
spectively. 
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IV. MEAN FIELD THEORY 



We focus on repulsive interactionis with [/ > in 
this paper. According to Landau theory, the condition 
(0, 0) = may be used to determine the critical point 
of the phase transition between the superfluid and the 
normal liquid Q . 

The key difference between the slave boson and slave 
fermion is their quantum statistics, which leads to the 
sign difference =F in equation (18). The difference ap- 
pears because we have approximated all Ai(r) by a real 
constant A. We discuss the zero temperature and the 
finite temperature cases separately. 



T/zt 

c 

1.5 . 




A. Zero Temperature 



In the zero temperature limit, Dickersheid et al j3| in- 
vestigated commensurate fillings and assume the num- 
ber of particles of each well to be fixed at some value a': 
n" — 5a^a' in mean field theory. This mean field assump- 
tion works for both kinds of slave particles. In terms of 
G~^(0, 0) = and the Green's function (|I7|I . it is easy 
to calculate the phase boundaries in the ^i-U plane [g: 
The Mott insulator phase is in the regimes where /2 lies 
between ft'^ 

ni - \[Ui2a' - 1) - 1] ± iy^C72-2C7(2a' + I) + l(.19) 

Here fl — [ijzt and U — U / zt. It is easy to check that 
for U > 0, ^ is positive. Eq. (fT^ reproduce results of 
previous mean field studies [9|. For the first Mott lobe, 
a' = I, the critical temperature is give by the zero of 
the square root in ((T^ . which gives tJc ~ 5.83. The 
quantum Monte Carlo calculation showed the critical Uc 
is somewhat smaller: C/c « 4 [l^. 



B. Finite Temperature 



FIG. 2: The critical temperature of the normal-superfluid 
phase transition as a function of the interaction U . The solid 
curve is for the ct = 0,1,2 slave fermions. The dash curve 
is for the a = 0,1,2 slave bosons. The down-triangles for 
the a — 0,1,2,3 slave fermions. The up-triangles for the 
a = 0,1,2,3 slave bosons. The data for the slave bosons 
are taken from The arrows indicate the position of the 
maxima of the critical temperature. The longer bar on the 
Tc-axis is the critical temperature (r^'^^'^'- = 1.18) for the ideal 
Bose gas in three-dimensional lattice. 
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We next turn our attention to the finite temperature 
behavior. For an site-independent parameter A, the re- 
laxed constraint associated with (|3Jl may be derived by 
taking a saddle point approximation for A, i.e., minimiz- 
ing the thermodynamic potential: dfl/dX — 0. The par- 
ticle conservation condition reads —dVl/dn = N . The 
corresponding saddle point equations are given by 



Nsl^{l - " ) - 7^ 2^ ^(k, iw„) = 0, 



iV, ^ an" + ^ ^ G(k, iLu^ 



9G-i(k,iL^„) 



k.n 



iV.(20) 



In these equations, all the possible multi-occupant states 
are included. 



FIG. 3: The a = 1 slave fermion number versus U/zt at 
the critical point of the phase transition. The solid curve is 
for three slave fermions and the dash line is for four slave 
fermions. The empty circles represent the unphysical results 
with > 1 at the "critical temperature". 



The mean field approximation implies the last terms in 
the above two equations may be neglected, which gives 



1, and 



E 



an 



N/L. 



(21) 



From the first equation, we see that all multi-occupancy 
of the different types slave particles are excluded in the 
mean field approximation |l4| . The difference between 
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the slave boson and slave fermion approaches is shown 
only on their different statistics. The critical point of 
the superfluid/normal liquid phase transition, in terms 
of G~^{0,0) = 0, is determined by 



/i — all 



(22) 



We now restrict to the commensurate state and take the 
total density to n = N/L — 1. To be able ito solve 
eqs. H21II and H22(l . one has to make a cut-off in a. For 
a cut off a A/ — 2 so that only a = 0,1,2 are allowed, 
these equations can be analytically solved. The resulting 
equation have been solved by Dickerscheid et al in the 
slave boson case They are also easily solvable for 
the slave fermion case. For a given U, fi = U/2 and 
= {U+3)/9 for both kinds of slave particles the critical 
point is determined by 



2 



U 

"2 



In 



In 



(C7-24)(C7 + 3) 
(C7-6)(C/ + 12) 

(£7 + 12)(i7 + 3) 
(C7-6)2 



(boson) (23) 



(fermion) (24) 



where Tc = Tc/zt. There is no analytical solution when 
ttM > 3, but the equations may be numerically solved. 
The results of four types of the slave bosons have been 
presented in Ref. We have solved the slave fermion 
equations for aM = 3. In Fig. 2, we plot the critical 
temperature as a function of U. Below Tc, it is superfluid 
phase and above Tc, it is normal liquid phase. Tc = is a 
triplet critical point of the superfluid, normal liquid and 
Mott insulator phases. 

We now make a comparison with the results in Fig. 2 
from the two approaches. The critical temperature for 
a given U from the slave boson approach is higher than 
that from the slave fermion one. Two further analyzes 
can show that the slave fermion result may be more ap- 
propriate. The bar located at T = 1.18 on the T-axis is 
the exact critical temperature for the ideal Bose gas on 
the three-dimensional lattice in the long wave length limit 
[l5|. It is obvious that the critical temperature at [/ = 
from the slave fermion approach is better then that from 
the slave boson approach. As mentioned in Q , there are 
two unsatisfactory features for slave bosons: one is that 
the zero-temperature critical interaction U is moved to 
Uc ~ 6 from solving (|21|) and (|22|l and is different from 
zero-temperature result Uc ~ 5.83 from given by (|19|l . 
Another is the existence of a local maximum in the po- 
sition of the phase boundary indicated by the arrows in 
Fig. 2. We see that both features are improved by the 
slave fermion approach. Although the three slave fermion 
result still gives J7c = 6 at T = 0, the four fermion result 
has moved Uc to about 5.85, which is very close that the 
well-known zero-temperature mean field result Uc ~ 5.83. 
Note that there are solutions of J7 > 6 in eq. H24() where 
the 'critical temperatures' are nonzero, which are shown 
by filled circles in Fig. 2. However, these solutions cor- 
respond to > 1 and n° < 0. Thus, these points are 



not physical. The similar situation appears for four slave 
fermions but these points move to J7 > 5.85 (shown by 
the empty circles in Fig. 2). To verify these 'critical tem- 
peratures' are not physical, we depict at the critical 
point as a function of U in Fig. 3. We see that > 1 
after U > 5.85. All of these evidences show the critical 
point indeed moves to around Uc ~ 5.85 at T = 0, which 
is in good agreement with the zero temperature mean 
field result, Uc = 5.83. 

The second feature in the slave boson result is the lo- 
cal maximum in the U-Tc curves (see the arrows in Fig. 
2). As we see soon, this maximum comes from the finite 
aM approximation getting worse as J7 — > 0. This fea- 
ture is greatly improved in the slave fermion U-Tc curves. 
The positions of the maximum of the critical tempera- 
ture for slave fermions are in about U ^ 1.2 for three 
fermions and near U = 0.2 for four fermions, which are 
much closer to zero than 2.15 and 1.8 in the slave boson 
curves. Furthermore, AT = f^^"" - fc{U = 0) ^ 0.25 
and AT/T™'^'' - 15% for slave bosons while AT ~ 0.02 
and AT/T™'"'" ^ 2% for slave fermions. The local maxi- 
mum appears because the cut-off a < 3 is not appropriate 
as [/ ^ 0. 

Both these two quantitative improvements from the 
slave fermion approach over the slave boson approxima- 
tion come from the exclusion of the same types of the 
slave fermions, due to the Pauli principle. From H23|l and 
(Pl|l . we can calculate the derivative It is seen that 
the boson statistics of the slave boson sharpens the slope 
the Tc-U curve in small U. 

We next investigate if the finite aM approximation 
is good or not. For this purpose, we plot n" {aM — 3) 
versus U at the critical temperature (Fig. 4). It is seen 
that the occupancy of the 0th, 1st and 2nd types of 
the slave fermion is of the order 10^^ from U ~ 1 to 
4. However, n'^ decreases quickly as U increases. For 
[/ = 1, n3 « 0.05 but « 0.009,5 x lO"'' and 2 x 10"^ 
for [/ = 2, 3 and 4, respectively. This means that for a 
large enough U, the a^i — 3 cut-off is a good approx- 
imation. In the regime of small U, a larger aAf(> 3) 
is required if we would like to have a quantitatively 
reliable result. For aM = 3, we see that, from Fig. 2, 
Tc{U — 0, aM — 3) ~ 0.98. It may be expected that 
as aM increases, Tc{U = 0) should be close to, e.g, 
1.18 in three dimensions, the ideal Bose gas critical 
temperature. The approximation getting worse for small 
U means the contribution from large a(> 3) can not be 
neglected. The maximum of the critical temperature in 
Fig. 2 comes from neglecting these degrees of freedom 
corresponding to large a. Taking a larger aM, it may 
be anticipated that the maximum of Tc may disappear 
as Tc{U — 0) tends to 1.18. To reveal the quantitative 
behavior of the system for small U more precisely, we 
have to work at a larger aM ■ The numerical work is still 
in progress which will be present elsewhere. 
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Incommensurate State 



0.2- 




1.0 1.5 2.0 2.5 3.0 3.5 4.0 
U/Zt 



It was also known that in the incommensurate state, 
there is no Mott insulator phase for any value of U. To 
confirm this point, we calculate two incommensurate fill- 
ings with n = 0.9 and 0.75. As Fig. 5 shown, only 
a normal-superfluid transition in the phase diagram is 
found and there is no Mott insulator phase. 



CONCLUSION 



FIG. 4: The slave fermion occupant number versus U at 
the critical temperature (oa/ = 3). 




0.4- 



2 4 6 8 10 

U/zt 



FIG. 5: The phase diagram for the incommensurate fillings. 
The filled squares are for n — 0.9 and the empty circles are 
for n = 0.75. 



In summary, we have made a comparative study of the 
finite temperature phase diagrams with the slave boson 
and the slave fermion mean field theory for the Bose Hub- 
bard model. We found that both slave particle mean field 
theories are qualitatively the same but the slave fermion 
approach is quantitatively more accurate. Many other re- 
sults obtained in Ref. Q by the slave boson approach are 
valuable and may be improved by the slave fermion ap- 
proach by replacing the bosonic occupation number with 
the the fermionic occupation number. This approach 
can be generalized to the two-component Bose Hub- 
bard model with 6,^^ = Ea^^a^ V^T + lcl^+i,„^CaT,ai 

and = Ea^.a^ V'^H^^T.^i-nCaT.^i' "^^il^ for 

Bose-Fermi mixture Hubbard model, a slave boson- 
composite fermion mixture approach, with b] = 

Ea Va + l(cl,a+lC»a+^>+l&»a) and // = Ea 4,aho^ OT 

~ Ea ^1 a'^i,ai "^iH bc required. Here cj ^ may be thought 
as a composite fermion (1^. These will be worked out 
elsewhere. 

This work was supported in part by Chinese National 
Natural Sciences Foundation and NSF of the US. 
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